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Article history: An electroencephalogram (EEG) is a recording of various frequencies of 
. electrical activity in the brain. EEG signal is very useful for diagnosis of 
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paper, an efficient algorithm is proposed to process the raw EEG signal to 


combat the noise. To obtain noiseless EEG data, the likelihood test ratio is 
Keywords: applied to interference computation block. The likelihood ratio test converts 
EEG data signal into segmented data with nearly constant noise 


San characteristics. This will aid in detecting the noise present in a tiny segment 
enolsing which ensures proper signal denoising. The processed signal is compared 
Electroencephalogr. —_ with the database of noiseless EEG of the same person using principal 
Interference calculation component analysis (PCA) classifier. The proposed algorithm is 99.01% 
Principle component analysis efficient to identify and combat noise in the EEG signal. 
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1. INTRODUCTION 

Recently after stroke, epilepsy is the second most severe neurological disorder viewed in human. 
Nearly 50 million people, or around 1% of the global population, are affected. One out of every ten people 
will experience a seizure at some point in their life, but the majority of them will not have epilepsy because 
the root cause of the convulsions is not related to the brain. At any point in their life span, one in every 50 
people will develop epilepsy. Everyday, about 75 new cases of epilepsy are discovered. Flashing lights can 
affect just around 3-5% of epilepsy sufferers (photosensitive epilepsy). It’s a common myth that all epileptics 
are affected by rapid visual stimuli; however, this is not the case. This disease affects just about two out of 
every 10,000 people in the general population [1]. It is estimated that 1.5 to 2 million people in the United 
States alone suffer from active epilepsy at this time [2]. Many people with epilepsy today have lives that have 
been significantly changed by medical science. Indeed, many people with epilepsy have gone on to have 
good careers or become very well-known. Alfred Nobel, the inventor of dynamite and the founder of the 
Nobel Prize, suffered from epilepsy. Hans Berger, a German psychiatrist, invented the electroencephalogram 
(EEG) test in 1924 and demonstrated its functional and diagnostic utility. Sheng Li and Hanxin Feng had 
proposed EEG signal classification method based on feature priority analysis and convolutional neural 
network (CNN). The importance of detection electrodes is sorted by random forest algorithm, and the higher 
priority electrodes are screened out. Hence, they concluded that proposed method is effectively realizes 
feature extraction and classification of EEG signals [3]. Special sensors are strategically positioned [4] 
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around the head and linked to a computer that tracks electrical impulses on a screen or on paper. Trained 
neurologists have examined the various frequencies in the EEG and identify patterns that provide details 
about the epileptic condition [5]. Alpha (7.5-13) Hz, beta (13-30) Hz, delta (0-3.5) Hz, and theta (3.5-7.5) Hz 
are the four basic brain waves commonly used to characterize raw EEG signals. These bands reflect the 
brain’s most active operation [6]. 

In this paper, effective noise minimization technique for real time EEG signal is proposed which is 
based on detecting the interferences using novel mathematical models. The interference calculation is used 
with likelihood test ratio detector to obtain the noise data small segmented EEG signal where noise variations 
remain constant. This helps in proper filtering supported by our obtained results through principal component 
analysis (PCA) classifier. 


2. LITERATURE SURVEY 

All Ictal states apply to the various stages of an epileptic seizure. In the broadest context, these 
states reflect the various phases of an epileptic seizure [7]. The following are the various stages of epilepsy. 
i) Interictal: a natural resting state without seizure disorder. 11) Preictal: a period of time preceding a seizure 
that does not refer to the brain’s normal state [8]. ii) Ictal state of seizure termination: The duration of the 
seizure’s activation phase [9]. iii) Postictal condition: The early strategies in the field were aimed at 
compressing data and, more importantly, PA-1-highlighting occurrences for subsequent evaluation by a 
neurologist. These semiautomated detection methods can expedite the review process, but also postpone 
intervention. Several papers [10]—[13] describe commonly used compression methods as compressed spectral 
array, density spectral array, spectrogram, and non-linear energy operator. The most extensively used 
strategy in the past was to evaluate the signal’s frequency content using a variation of the fourier transform. 
There is a significant trade-off between time and frequency resolutions when employing this strategy. The 
time resolution is good when evaluating the spectrum in short windows, but the frequency resolution is 
inadequate. The frequency resolution improves with a longer window, but the information becomes less 
concentrated in time. The wavelet transform is a more contemporary method. Wavelet decomposition allows 
you to depict the qualities of a signal on different scales. This way, both the time and the money are saved. 
Expert systems, decision trees, clustering algorithms, self-organizing maps, and a variety of artificial neural 
network topologies are examples of machine-learning methods that have been used in this subject. While 
advanced machine learning approaches can improve the algorithm’s performance, they can also be difficult 
to interpret for the end-user. The easier it is for specialists to learn to trust and understand a system, the 
simpler it is. Here are some of the most extensively utilised systems. Since the 1970s, Gotman has worked in 
the field of EEG monitoring and seizure detection. Stellate is the company that distributes his algorithms. The 
first algorithms examined the properties of EEG waves after decomposing them into constituent waves 
[14], [15]. A new module [16] was introduced to rule out typical reasons of false positive detections. Seizure 
detection software is also available from persyst development corporation. Reveal Rosetta, the algorithm, is 
also expected to have ICU applications. The algorithm’s structure has mostly remained unknown. In [17] 
causal network elicitation technique (CNET) is a non-commercially distributed seizure detection programme. 
It describes cepstral characteristics [18]. 


3. PROPOSED MODEL 

The flow diagram of proposed method is shown in the Figure | which consists of proposed 
de-noising method, PCA classifier and comparator blocks. The different types of noises are added by the 
noise generator algorithm present in matrix laboratory (MATLAB) tool. The interference calculation block is 
used to calculate the noise and interferences present in the signal using proposed mathematical modelling. 
The noise minimization block is used to detect noises and perform smoothing of the signal through proposed 
mathematical model. Now using the separate PCA signal classification [19]—[21] algorithms the actual EEG 
signal (noiseless) and smoothen EEG signal are classified into different state’s which are compared using 
eucledean distance to check the effectiveness of the de-noising algorithm. 


3.1. Read EEG signal 

The EEG signals are read through the “Read EEG Signal” block which contains certain artifacts. 
Artifacts can originate from various sources such as the subject, equipment, or the environment and consist of 
ocular artifacts, such as eye blinks; movement of the EEG sensors; and electromyogenic artifacts, caused by 
muscle movement. The standard database [22] are used for the calculation. 
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Figure 1. Proposed noise reduction method with performance analysis 


3.2. Interference calculations 

Dynamic interference calculation method is used to detect the artifacts. In this method, interferences 
are analysed in learned subspace followed by a novel technique generated from modified robust subspace 
detection [1] method. In this method, first the unknown interferences are calculated and then using those 
interferences the actual interferences are calculated. Let us consider the EEG signal consists of EEG data 
with some unknown noise which is expressed mathematically as 


x=rnt+ntn (1) 


where, r; is the actual EEG data, r, is the unknown interference, and 1 is the noise. The (1) is for the space 
containing the entire EEG signal. To get accurate interference, subspace calculation method is used. Let us 
consider S is the subspace which is used to estimate the data. As a result, K vector and M-dimensional 
subspace is used to define the total subspace S which modifies the (1) as 


x= SO+UP+n (2) 


where, @ is the unknown gain of the noise, U is noise amplitude, and @ is phase of the noise. Through 
likely-hood ratio test [1] of (2), the equation for log yields become 


Ae) =(L|x-sd],) + (Ale -all,) (3) 


The use of likely hood ratio test for the (2) results in segmentation of entire EEG data signal into very small 
segment where the noise characteristics are almost constant. This will helps to track the proper noises present 
in small segment which in turn helps to perform proper de noising of the signal. Since the feature vector x to 
be in a space with components in the direction of and orthogonal to the signal subspace S, the unknown gain 
of the signal is 
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sist 
Ss "sts (4) 
By considering the ratio of the portion of x in the subspace of S, the interference equation is then 
XT(Ps—Pn)X 
A(x) = See (5) 
_ XT Ps.X 
Mx) = XT Py.X (6) 


where, Pn is projection matrix onto the part of the measurement space orthogonal to the signal. 


3.3. Noise minimization of EEG signal 

To minimize noises, it is essential to calculate the noises present in the input signal. In this case, the 
interference equation explained in the “Interference Calculation” block are used to learn the interference 
characteristics with respect to the calculated inferences. The “Adaptive Matched Subspace detector” [2] is 
then used to minimizes these noises which is written as 


XT (Ps—Pp)X 
XT(Py—Pp)X 


A(x) = (7) 
3.4. Classification of noise 
The general PCA [3] technique is used to classify the signals. 


a) (8) 


Oo 


where, x symbols and pu mean value. The covariance matrix equation is then 


06 6 ———— (9) 


OD Tee 
The final data is then 
final data = feature vector * transpose(scaled(data)) (10) 


3.5. Compare 

The eucledian distance is used to compare the features. Now by using different threshold values, the 
ROC is calculated. The different threshold values are used to calculate the corresponding False acceptance 
ratio (FAR) and total success rate (TSR) values [23] using separately noisy and noiseless signals from the 
(10). Similarly the accuacy calculation equation is given as 


Accurate —Denoised 


Accuarcy = ( ) x 100 (11) 


Accurate 


4. SIMULATION RESULT AND DISCUSSIONS 

The proposed algorithm is implemented on MATLAB tool where the entire algorithm is coded using 
standard MATLAB programming techniques. The standard EEG database [22] are used to check the 
performance of the algorithm. To calculate the detection accuracy, the entire databases are used. But for 
simplicity, only the snapshots of three persons are shown. 


4.1. Case I: person 1 

The resulting receiver operating characteristic (ROC) graph is shown in the Figure 2, it can be seen 
that the both graph of our proposed method for smoothened version of noisy EEG data and the existing 
noiseless present in the EEG database for the person | is almost similar and is overlapped which proves the 
effectiveness of the algorithm. The analysis of the detection of interference with changing magnitude is 
plotted on fixed threshold at fixed noise level is shown in the Figure 3 for existing technique [24] mentioned 
in blue color and proposed method mentioned in red color. It can be seen that True positive rate is same for 
both technique with the proposed technique is nearly generating same curve almost similar to noiseless 
signal. 
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Figure 2. ROC curve of actual and detected interferences of database image 1 
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Figure 3. Interference magnitude of actual and detected interferences for database image 1 


The analysis of the detection of interference with changing magnitude is plotted on fixed threshold 
at fixed noise level is shown in the Figure 3 for existing technique [24] mentioned in blue color and proposed 
method mentioned in red color. It can be seen that true positive rate is same for both technique with the 
proposed technique is nearly generating same curve almost similar to noiseless signal. The false positive rate 
with interference magnitude is plotted in Figure 4 and it is observed that the proposed method shown in red 
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the algorithm and it is visible that the dynamic tracking mechanism can track interference more efficient than 
existing. Since small duration characteristics are defined properly. 
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Figure 4. Detection of actual and detected interferences at fixed threshold and fixed noise database image | 


4.2. Case II: person 2 

Similar to Figure 2, the ROC curve of different for different random person (person 2) present in 
standard database [22] is considered to prove the efficiency as shown in the Figure 5 which shows that the 
both graph of proposed technique for our smoothened version of noisy EEG and existing noiseless in 
database is almost similar and overlapped which proves the effectiveness of the algorithm. The analysis of 
the detection of interference with changing magnitude is plotted on fixed threshold at fixed noise level is 
shown in the Figure 6 for existing technique [24] mentioned in blue color and proposed method mentioned in 
red color. It can be seen that true positive rate is same for both technique with the proposed technique is 
nearly generating same curve almost similar to noiseless signal. The false positive rate with interference 
magnitude is plotted in Figure 7 and it is observed that the proposed method shown in Red color is minimum 
compared to existing [24] technique shown in blue color which shows the effectiveness of the algorithm and 
it is visible that the dynamic tracking mechanism can track interference more efficient than existing. Since 
small duration characteristics are defined properly. 
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Figure 5. ROC curve of actual and detected interferences of database image 2 
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Figure 6. Interference magnitude of actual and detected interferences for database image 2 
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Figure 7. Detection of actual and detected interferences at fixed threshold and fixed noise database image 2 


4.3. Case III: person 3 

Similar to Figure 2, the ROC curve of different for different random person present in standard 
database [22] is considered to prove the efficiency as shown in the Figure 8 which shows that the both graph 
of proposed technique for our smoothened version of noisy EEG and existing noiseless in database image 3 is 
almost similar and overlapped which proves the effectiveness of the algorithm. The analysis of the detection 
of interference with changing magnitude is plotted on fixed threshold at fixed noise level is shown in the 
Figure 9 for technique [24] mentioned in blue color and proposed method mentioned in red color. It can be 
seen the true positive rate is same for both technique with the proposed technique is nearly generating same 
curve almost similar to noiseless signal. The false positive rate with interference magnitude is plotted in 
Figure 10 and it is observed that the proposed method shown in red color is minimum compared to existing 
[24] technique shown in blue color which shows the effectiveness of the algorithm and it is visible that the 
dynamic tracking mechanism can track interference more efficient than existing. Since small duration 
characteristics are defined properly. 
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Figure 10. Detection of actual and detected interferences at fixed threshold and fixed noise database image 3 
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5. COMPARISONS WITH EXISTING TECHNIQUES 

The accuracy is compared in the below Table 1, which shows that the proposed technique is better 
in terms of accuracy with respect to existing technique presented by Balamareeswaran and Ebenezer [25]. 
The discrete wavelet transform (DWT) is used to denoise input EEG signal and the weighted support vector 
machine (SVM) technique is used for classification. On other hand the use of likekihood ratio test in our 
proposed method converts EEG data signal into segmented data of constant noise charactristics for proper 
denoising . Further, the denoised signal is used to characterize the disease depending upon variable threshold 
using PCA classifier to achieve higher accuracy rate. 


Table 1. Comparison of existing technique with proposed technique 


Authors Technique Accuracy 
Balamareeswaran and Ebenezer [25] DWT 97.66% 
Proposed Statistical Mathematical Model 99.01% 


6. CONCLUSION 

An efficient de-noising algorithm for EEG signal is proposed in this paper where the noise present 
in the raw EEG signal is detected through novel mathematical modelling and depending upon the noise 
characteristics the de-noising has been performed. The only limitation is if the noise is fully random then the 
tracking may not be effective. Since the algorithm segments the entire EEG signal in to smaller segments 
depending on similar or same kind of noise characteristics. The de-noised signal is then used to characterize 
the disease depending upon variable threshold using general PCA classifier. The above graph is compared 
with the graph generated by corresponding noise free EEG signal using same technique. This shows that the 
proposed algorithm can pre-process the signal effectively due to effective novel mathematical tracking of 
noises. 


REFERENCES 

[1] B. F. Shneker and N. B. Fountain, “Epilepsy,” Disease-a-Month, vol. 49, no. 7, pp. 426-478, Jul. 2003, doi: 10.1016/S0011- 
5029(03)00065-8. 

[2] I.E. Leppik, Epilepsy: a guide to balancing your life. New York: Demos Medical Publishing, LLC, 2007. 

[3] S. Li and H. Feng, “EEG signal classification method based on feature priority analysis and CNN,” in 2019 International 
Conference on Communications, Information System and Computer Engineering (CISCE), Jul. 2019, pp. 403-406, doi: 
10.1109/CISCE.2019.00095. 

[4] L. M. Hively and V. A. Protopopescu, “Channel-consistent forewarning of epileptic events from scalp EEG,” IEEE Transactions 
on Biomedical Engineering, vol. 50, no. 5, pp. 584-593, May 2003, doi: 10.1109/TBME.2003.8 10693. 

[5] J. Aicardi, M. A. Dichter, J. Engel, T. A. Pedley, E. Perucca, and M. R. Trimble, Epilepsy: a comprehensive textbook. Wolters 
Kluwer/Lippincott Williams and Wilkins, 2007. 

[6] S. Sanei and J. A. Chambers, EEG signal processing. West Sussex, England: John Wiley and Sons Ltd, 2007. 

[7] A. Pitkéanen, P. S. A., and S. L. Moshé, Models of seizures and epilepsy. Elsevier, 2006. 

[8] J. Gotman, “Automatic detection of epileptic seizures,” in Handbook of Clinical Neurophysiology, Elsevier, 2003, pp. 155-165. 

[9] J. Gotman, “Automatic detection of seizures and spikes,” Journal of Clinical Neurophysiology, vol. 16, no. 2, pp. 130-140, Mar. 
1999, doi: 10.1097/0000469 1 -199903000-00005. 

[10] M. L. Scheuer, “Continuous EEG monitoring in the intensive care unit,” Epilepsia, vol. 43, pp. 114-127, Jun. 2002, doi: 

10.1046/j.1528-1157.43.s.3.7.x. 

[11] R. Agarwal and J. Gotman, “Long-term EEG compression for intensive-care settings,” IEEE Engineering in Medicine and 

Biology Magazine, vol. 20, no. 5, pp. 23-29, 2001, doi: 10.1109/51.956816. 

[12] M. L. Scheuer and S. B. Wilson, “Data analysis for continuous EEG monitoring in the ICU: seeing the forest and the trees,” 

Journal of clinical neurophysiology: official publication of the American Electroencephalographic Society, vol. 21, no. 5, pp. 

353-378, 2004. 

[13] A. K. Shah, R. Agarwal, J. R. Carhuapoma, and J. A. Loeb, “Compressed EEG pattern analysis for critically Ill neurological- 

neurosurgical patients,” Neurocritical Care, vol. 5, no. 2, pp. 124-133, 2006, doi: 10.1385/NCC:5:2:124. 

[14] J. Gotman, “Automatic recognition of epileptic seizures in the EEG,” Electroencephalography and Clinical Neurophysiology, vol. 

54, no. 5, pp. 530-540, Nov. 1982, doi: 10.1016/0013-4694(82)90038-4. 

[15] J. Gotman, “Automatic seizure detection: improvements and evaluation,” Electroencephalography and Clinical Neurophysiology, 

vol. 76, no. 4, pp. 317-324, Oct. 1990, doi: 10.1016/0013-4694(90)90032-F. 

[16] H. Qu and J. Gotman, “Improvement in seizure detection performance by automatic adaptation to the EEG of each patient,” 

Electroencephalography and Clinical Neurophysiology, vol. 86, no. 2, pp. 79-87, Feb. 1993, doi: 10.1016/0013-4694(93)90079- 

B. 

[17] S. B. Wilson, M. L. Scheuer, R. G. Emerson, and A. J. Gabor, “Seizure detection: evaluation of the Reveal algorithm,” Clinical 

Neurophysiology, vol. 115, no. 10, pp. 2280-2291, Oct. 2004, doi: 10.1016/j.clinph.2004.05.018. 

[18] A. J. Gabor, R. R. Leach, and F. U. Dowla, “Automated seizure detection using a self-organizing neural network,” 

Electroencephalography and Clinical Neurophysiology, vol. 99, no. 3, pp. 257-266, Sep. 1996, doi: 10.1016/0013- 

4694(96)96001-0. 

[19] F. Tahernezhad-Javazm, V. Azimirad, and M. Shoaran, “A review and experimental study on the application of classifiers and 
evolutionary algorithms in EEG-based brain—machine interface systems,” Journal of Neural Engineering, vol. 15, no. 2, Art. no. 
021007, Apr. 2018, doi: 10.1088/1741-2552/aa8063. 


Efficient de-noising technique for electroencephalogram signal processing (Virupaxi Dalal) 


o ISSN: 2252-8938 


[20] Y. Lu, H. Jiang, and W. Liu, “A review of EEG signal classifier based on deep learning,” in Proceedings of Information Science 
and Cloud Computing-PoS(ISCC 2017), Feb. 2018, Art. no. 60, doi: 10.22323/1.300.0060. 

[21] A. S. R. and H. C. Nagaraj, “Classification of EEG signal using EACA based approach at SSVEP-BCI,” IAES International 
Journal of Artificial Intelligence (IJ-Al), vol. 10, no. 3, pp. 717—726, Sep. 2021, doi: 10.1159 1/ijai.v10.i3.pp7 17-726. 

[22] “EEG database,” Available: https://physionet.org/content/eegmat/1.0.0/, June 2021. 

[23] S. S. Bhairannawar, S. Sarkar, K. B. Raja, and K. R. Venugopal, “Implementation of fingerprint based biometric system using 
optimized 5/3 DWT architecture and modified CORDIC based FFT,” Circuits, Systems, and Signal Processing, vol. 37, no. 1, pp. 
342-366, Jan. 2018, doi: 10.1007/s00034-017-0555-0. 

[24] G. Li and J. J. Jung, “Maximum marginal approach on EEG signal preprocessing for emotion detection,” Applied Sciences, vol. 
10, no. 21, Art. no. 7677, Oct. 2020, doi: 10.3390/app10217677. 

[25] M. Balamareeswaran and D. Ebenezer, “Denoising of EEG signals using discrete wavelet transform based scalar quantization,” 
Biomedical and Pharmacology Journal, vol. 8, no. 1, pp. 399-406, Jun. 2015, doi: 10.13005/bpj/627. 

BIOGRAPHIES OF AUTHORS 


Virupaxi B. Dalal © i P is a research scholar at Department of Electronics and 
Communication Engineering, SDMCET, Dharwad, Karnataka, India. Author is also working as 
Assistant Professor in Department of Electronics and Communication at Jain College of 
Engineering and Research, Belagavi, Karnataka, India. Author’s area of intrest are Signals 
Processing and Image Processing. He can be contacted at email: virupaxi.dalal @ gmail.com. 


Dr. Satish S Bhairannawar © —4 P holds the position of Dean of Industry Institute 
Interface (III), Chairperson, startup incubations and Professor of Electronics and Communication 
Engineering, SDMCET, Dharwad. He has been involved in research activities funded by reputed 
state and central govt. organizations and has published over 22 research publications in SCIE. He 
can be contacted at email: satishbhairannawar @ gmail.com. 


Int J Artif Intell, Vol. 11, No. 2, June 2022: 603-612 


